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Abstract. Long-slit NTT-I-EMMI echellograms of NGC 6818 (the Little Gem) at nine equally spaced position 
angles, reduced according to the 3-D methodology introduced by Sabbadin et al. (2000a,b), allowed us to derive: 
the expansion law, the diagnostics and ionic radial profiles, the distance and the central star parameters, the 
nebular photo-ionization model, the 3-D reconstruction in He II, [O III] and [N II], the multicolor projection 
and a series of movies. The Little Gem results to be a young (3500 years), optically thin (quasi-thin in some 
directions) double shell (Mion—0.13 M©) at a distance of 1.7 kpc, seen almost equatorial on: a tenuous and patchy 
spherical envelope (r~0.090 pc) encircles a dense and inhomogeneous tri-axial ellipsoid (a/2~0.077 pc, a/b~1.25, 
b/c~1.15) characterized by a hole along the major axis and a pair of equatorial, thick moustaches. NGC 6818 
is at the start of the recombination phase following the luminosity decline of the 0.625 Mq central star, which 
has recently exhausted the hydrogen shell nuclear burning and is rapidly moving toward the white dwarf domain 
(log T* ~5.22 K; log Lt/L© ~3.1). The nebula is destined to become thicker and thicker, with an increasing 
fraction of neutral, dusty gas in the outermost layers. Only over some hundreds of years the plasma rarefaction 
due to the expansion will prevail against the slower and slower stellar decline, leading to a gradual re-growing 
of the ionization front. The exciting star of NGC 6818 (my —17.06) is a visual binary: a faint, red companion 
(my ~17.73) appears at 0.09 arcsec in PA=190°, corresponding to a separation >150 AU and to an orbital period 
>1500 years. 
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1. Introduction 

It is widely accepted that an “aged” star of low to medium 
mass (1.0 M0<Mms<8.O Mq) which is evolving toward 
the white dwarf region of the H-R diagram, first expels 
and then ionizes the surface layers, thus generating the 
Planetary Nebula (PN) phenomenology (Paczynski 1970, 
Aller 1984, Pottasch 1984, Osterbrock 1989). 

Thanks to the advent of sophisticate photo-ionization 
codes (Harrington 1989, Pequignot 1997, Ferland et al. 


* Based on observations made with ESO Telescopes at the La 
Silla Observatories, under programme ID 65.1-0524, and on ob¬ 
servations made with the NASA/ESA Hubble Space Telescope, 
obtained from the data archive at the Space Telescope Institute 
(observing programs GO 7501 and GO 8773; P.I. Arsen 
Hajian). STScI is operated by the association of Universities 
for Research in Astronomy, Inc. under the NASA contract NAS 
5-26555. We have applied the photo-ionization code CLOUDY, 
developed at the Institute of Astronomy of the Cambridge 
University. 
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1998), we now know in detail the physical effects produced 
by any UV flux on any gas distribution and composition. 
Conversely, the situation is quite disappointing for a true 
nebula, due to projection limitations: on the one hand the 
HST imagery has strongly enhanced the ground-based ev¬ 
idences of the PNe complexity, on the other hand the ob¬ 
servational data are still interpreted in terms of approxi¬ 
mate structures and unrealistic assumptions for the physi¬ 
cal conditions, like Te=constant and iVe=constant all over 
the object (more comments are in Aller 1984, 1990, 1994). 

In order to overcome the wide gap between theory 
and practice, the apparent, bi-dimensional nebular image 
should be de-projected, and the accurate spatial distribu¬ 
tion of the gas recovered. 

To this end we have developed an original procedure 
based on high dispersion spectra: the PN being an ex¬ 
tended and expanding plasma, the position, thickness and 
density of each elementary volume can be, in principle, ob¬ 
tained from the radial velocity, width and flux of the corre¬ 
sponding emission. We first apply a tomographic analysis, 
which reconstructs the ionic distribution in the nebular 
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This figure is available in JPG format 


Fig. 1. WFPC2 archive images of NGC 6818 in [O III] (left) and [N II] (right). The binary nature of the central star 
is shown in the enlargement at the bottom. North is up and East to the left. 


slices covered by the spectrograph slit, and then assemble 
all the tomographic maps by means of a 3-D rendering pro¬ 
cedure for studying the morphology, physical conditions, 
ionization, spatial structure and evolutionary status. 

The rationale of the method and the earliest, rough re¬ 
sults based on plate echellograms go back to Sabbadin et 
al. (1985, 1987). More recent, “quantitative” observations 
(i.e. using a linear detector) concern NGC 40 and NGC 
1501, both objects covered at moderate spectral resolu¬ 
tion, R~ 20000-25000, with the Echelle -I- 1.82m telescope 
of Padua Observatory at Asiago, Cima Ekar (Sabbadin et 
al. 2000a, b, and Ragazzoni et al. 2001; hereafter Papers 
I to III, respectively). 

At the same time we have carried out a sur¬ 
vey of two dozen PNe and proto-PNe in both hemi¬ 
spheres with ESO NTT-I-EMMI (spectral range AA3900- 
7900 A, R=60000, spatial resolution~1.0 arcsec) and 
the Telescopic Nazionale Galileo (TNG)-I-SARG (spectral 
range AA4600-8000 A, R= 115000, spatial resolution~0.7 
arcsec). The observed sample covers a variety of morpholo¬ 
gies, kinematics and evolutionary phases, including NGC 
2392 (the Eskimo nebula), NGC 3132 (the Eight-Burst 
n.), NGC 3242 (the Ghost of Jupiter n.), NGC 6210 (the 
Turtle n.), NGC 6543 (the Cat’s Eye n.), NGC 6751 (the 
Glowing Eye n.), NGC 6826 (the Blinking n.), NGC 7009 
(the Saturn n.), NGC 7662 (the Blue Snowball n.), IC 
418 (the Spirograph n.). He 2-47 (the Starfish n.), MyCn 
18 (the Hourglass n.), MZ 3 (the Ant n.), and as many 
un-dubbed, but equally exciting targets. 

The 3-D ionization structure of NGC 6565 has been 
presented by Turatto et al. (2002, Paper IV); here we dis¬ 
cuss the case of NGC 6818 (PNG 025.8-17.9, Acker et al. 
1992). 

The paper is structured as follows: Sect. 2 introduces 
the nebula. Sect. 3 presents the observational material and 
the reduction procedure. Sect. 4 is dedicated to the gas 
kinematics, Sect. 5 concerns the radial profile of the phys¬ 


ical conditions (electron temperature and electron density 
from forbidden line ratios), in Sect. 6 the ionization struc¬ 
ture and the overall gas distribution are discussed, in Sect. 

7 we derive the nebular distance, mass and age, in Sect. 

8 the central star parameters are given. Sect. 9 contains 
the application of the photo-ionization model (CLOUDY), 
Sect. 10 describes the 3-D structure of the nebula in dif¬ 
ferent ions, and Sect. 11 presents a short discussion and 
the conclusions. 


2. The nebula 


The HST/WFPC2 appearance of the high excitation PN 
NGC 6818 (sometimes called the Little Gem) is shown in 
Fig. 1: in [O HI] “a roughly spherical outer envelope as 
well as a brighter vase-shaped interior bubble. There is a 
prominent orifice to the North and a smaller one to the 
South, along the major axis, probably caused by a blow¬ 
out from a fast wind” (Rubin et al. 1998). A few dark 
filaments and knots are also present. 

The [N H] emission mainly occurs in two irregular 
equatorial “moustaches” and in a multitude of condensa¬ 
tions (in some cases presenting a sort of radially arranged 
cometary tail), which are brighter in the southern part of 
the nebula. 


The striking HST multicolor reproduction by Arsen 


Hajian & Yervant Terzian (http://ad.usno.navy.mil/pne/ 


gallery.html) highlights the composite envelope structure 
and the large stratification of the radiation within NGC 
6818. 


Rubin et at. (1998) reported the presence of two faint 
stars, roughly 2-4 arcsec N and NE of the central star, 
which could be physically associated to it. A careful anal¬ 
ysis of the WFPC2 frames allowed us to unveil the true 
binary nature of the central star: the red, faint companion 
is at a distance of 0.093(±0.004) arcsec in position angle 
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(PA)= 190(±2)° (see the enlargement at the bottom of 
Fig. 1). 

Following Weedman (1968) our nebula is a prolate 
spheroid (an ellipse rotated about the major axis) of mod¬ 
erate ellipticity (axial ratio 3:2) observed perpendicular to 
the major axis, whereas Sabbadin (1984) suggests a some¬ 
what optically thick prolate spheroid in a early evolution¬ 
ary phase, seen almost pole-on. According to Hyung et al. 
(1999), a definitive statement on the geometry must wait 
for more adequate kinematical data, due to the complexity 
of NGC 6818. 

To this end we have secured a series of long-slit, radi¬ 
ally arranged echellograms. They were analysed with the 
3-D methodology described in Papers I to IV. 

3. The observational material and the reduction 
technique 

Nine echellograms of NGC 6818 (exposure time 600s; spec¬ 
tral resolution 60000 with a slit 1.0 arcsec wide) have been 
obtained on July 29, 2000 at the ESO NTT in photomet¬ 
ric sky conditions and seeing between 0.7 and 1.0 arcsec. 
We have used a 40 arcsec long slit centered on the excit¬ 
ing star for all the selected PA, ranging from 10° to 170° 
with a constant step of 20°. Since we do not insert an in¬ 
terference filter (as normally done by the other observers 
to isolate a single order), each spectrum covers 80 echelle 
orders (A3967 A of [Ne III] to A7751 A of [Ar III]), and 
provides the spatio-kinematical structure of a good two 
dozen nebular emissions, representing all the main ionic 
species. 

The reduction method follows conceptually the stan¬ 
dard procedure, including bias, flat held, distortion cor¬ 
rection, wavelength and flux calibration, and is carefully 
described in Paper IV. 

Fig. 2 illustrates the detailed structure in [O I] (A6300 
A), [N II] (A6584 A), H I (A6563 A), [O III] (A5007 A) 
and He II (A4686 A) at the observed PA. 

NGC 6818 exhibits a complex ionization structure: 

- [O I] is only seen in the outermost regions, in the form of 

distinct condensations, which in some cases are sym¬ 
metrically arranged (the moustaches), otherwise single 
(like the cometary knot in PA=150°, southern edge); 

- the low ionization [N II] emission forms a double en¬ 

velope of inhomogeneous structure from PA=50° to 
PA=130°, and a single, elongated and distorted ring 
from PA=150° to PA=30°. The presence of FLIERS 
(fast, low ionization emitting regions), as introduced 
by Balick et al. (1993, see also Corradi et al. 1996), 
is not conhrmed here: although some isolated spikes 
present a faint, high velocity tail, the overall [N II] 
emission of NGC 6818 can be understood in terms of 
“normal” nebular regions simply shadowed by some in¬ 
ner and dense layer causing the ionization drop in the 
outer plasma; 

- the Ha line (central row in Fig. 2) is characterized by 

a very blurred appearance (due to a combination of 


thermal motions plus fine structure plus expansion ve¬ 
locity gradient), masking the detailed distribution of 
the ionized gas. The same effect is also present in He I 
and He H, but it is particularly damaging for Ha, the 
reference emission in both the radial electron density 
and ionization determinations (Paper IV); 

- A5007 A of [O HI] (a mean-high excitation ion), by far 

the strongest line in the optical region, highlights the 
double envelope structure of NGC 6818: the inner shell 
(better seen between PA=50° and PA=130°) is an ir¬ 
regular ellipsoid broken along the major axis (N-S di¬ 
rection) whose equatorial, denser regions are identified 
by the “moustaches”. The line-tilt between PA=90° 
and PA=110° suggests that we are misaligned with 
both the intermediate and the minor axes of the inner 
ellipsoid, and that the line of the nodes is at PA~60°. 
The outer shell is almost un-tilted, spherical, holed at 
North and South, and circumscribes the internal one; 

- the He H emission (bottom row in Fig. 2) marks the 

highest excitation nebular regions, mainly constituted 
by the inner shell. Note the blurred appearance of 
A4686 A, essentially due to the thirteen hne structure 
components. 

The first, qualitative picture of the spatial structure com¬ 
ing from the echellograms confirms the indications already 
obtained from the imaging (Sect. 2). 

4. The gas kinematics 

According to Papers I to IV, the expansion velocity 
(Vexp) of the ionized gas can be derived from the analysis 
of the “central star pixel line” (cspl) in the different ions. 
The cspl is parallel to the dispersion, selects the nebu¬ 
lar material projected at the apparent position of the star 
(whose motion is purely radial) and is the same in all the 
frames, the slit being radially arranged. Thus, in order to 
improve the S/N of the faintest emissions, the nine echel¬ 
lograms have been combined. 

The results are contained in the last column of Table 
1, where the ions are put in order of increasing ioniza¬ 
tion potential (IP). Typical errors are 1.5 km s“^ for the 
strongest forbidden emissions (like A6584 A of [N H] and 
A5007 A of [O HI]) to 3.0 km s“^ for the faintest ones 
(in particular: A6300 A of [O I], because of the knotty 
structure, and A7319.87 A of [O H], which is partially 
blended with A7318.79 A, also belonging to the 0+ red 
quartet). The corresponding uncertainties for the recom¬ 
bination lines are: 2.5 km s“^ for A6563 A (H I), and 2.0 
km s“^ for A5876 A (He I) and A4686 A (He H). 

Golumns 3 to 6 of Table 1 report the kinematical re¬ 
sults from the literature. In detail: 

- Wilson (1950) obtained a single Coude spectrum of 

NGC 6818 (without de-rotator) at spectral resolution 
R~30000; 

- Sabbadin (1984) observed the nebula at four PA (long- 

slit echellograms at R~ 15000); 
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Fig. 2. Detailed structure (in logarithmic scale) of some representative emissions at the nine observed PA of NGC 
6818. The blue-shifted gas is to the left. The slit orientation at each PA is indicated in the corresponding [O I] frame. 
The A6300.304 A night sky line has been removed from the [O I] spectral image. 

Table 1. Peak separation in the cspl of NGC 6818 


Ion 

IP (eV) 



2Vexp (km/s) 



[0 I] 

0.0 

Wilson 1950 

Sabbadin 1984 

Meatheringham et al. 1988 

Hyung et at. 1999 

this paper 

72: 

[S II] 

10.4 

- 

- 

- 

60.7 

71 

[O II] 

13.6 

60.2 

- 

59.1 

58.9 

69: 

H I 

13.6 

55.5 

52 

- 

42.6 

55 

[N II] 

14.5 

- 

64 

- 

61.8 

70 

[S III] 

23.4 

- 

- 

- 

50.8 

65: 

He I 

24.6 

- 

- 

- 

57.3 

66 

]Ar III] 

27.6 

- 

- 

- 

48.9 

64 

[0 III] 

35.1 

56.2 

54 

55.1 

51.3 

62 

[Ar IV] 

40.7 

- 

- 

- 

30.7 

48 

]Ne III] 

41.0 

58.0 

- 

- 

52.2 

63 

N III 

47.4 

- 

- 

- 

41.5 

- 

He II 

54.4 

42.4 

- 

41.9 

31.5 

45 

[Ar V] 

59.8 

- 

- 

- 

20.2 

36 

[Ne V] 

97.1 

32.6 

- 

- 

- 

- 


- Meatheringham et al. (1988, long-slit echellograms at 

unspecified position, R~26000) also measured the 
width at 10% maximum intensity of A5007 A, obtain¬ 
ing 2Fexp[0 III]= 87.6 km s“^ (following Dopita et al. 
1985, this corresponds to the largest expansion velocity 
of the gas); 

- Hyung et al. (1999, R~33000) studied four bright re¬ 

gions located along the apparent minor and major axes 
(4-5 arcsec East and West, and 9-10 arcsec North and 
South of the central star); their values represent lower 
limits to Rexp(various ions), due to projection efi'ects. 

In Table 1 the kinematical variety reported by the dif¬ 
ferent authors is symptomatic of the difficulties connected 


to the analysis of the high dispersion spectra, and, accord¬ 
ing to Hyung et al. (1999), stresses the fact that only the 
detailed coverage at adequate spatial and spectral resolu¬ 
tions can provide a reliable information on the structure 
of a chaotic object like NGC 6818. 

As cspl refers to the kinematical properties of the mat¬ 
ter projected at the apparent position of the star, so the 
“zero velocity pixel column” (zvpc), corresponding to the 
recession velocity of the whole nebula, gives the radial 
distribution of the ionized gas which is expanding per¬ 
pendicularly to the line of sight (Paper IV and references 
therein). 

The intensity peak separations, 2rzvpc, in the different 
emissions at the nine PA of NGC 6818 are presented in 



S. Benetti et al.: The Planetary Nebula NGC 6818 


5 


Table 2. Both Fig. 2 and Table 2 indicate complex zvpc 
profiles, often multi-peaked (inner and outer shell), fast 
changing in direction, and variable (at a given PA) from 
ion to ion. The observed line structure depends on the 
radial matter distribution, the ionization and the intensity 
of the emission (the same occurs for the cspl assuming 
Pexpocr, as normally observed in PNe). 

To be noticed: the “relative” spectral and spatial 
resolutions of our echellograms, as introduced in Paper 
III, differ by a factor of two. The former is given by 
RR=Pexp/AV ~ 6, AV being the spectral resolution, 
and the latter by SS=r/Ar ~ 12 (r=apparent radius, 
Ar=seeing). In practice this means that the spatial infor¬ 
mation of NGC 6818 is twice as detailed as the kinematical 
one. 

In order to assemble the kinematical and the spatial 
results (corresponding to a single radial direction and to 
nine tangential directions, respectively) we must identify 
the PA at which r^vpa—rcspi- The solution comes from the 
qualitative picture resulting from both the imaging (Sect. 
2) and the spectra (Sects. 3 and 4). Let’s consider for the 
main component of NGC 6818 the most general spatial 
structure, i. e. a tri-axial ellipsoid. The major axis, pro¬ 
jected in PAc:ilO°, is almost perpendicular to the line of 
sight, since the emissions, although inhomogeneous and 
distorted by irregular motions, appear un-tilted. The in¬ 
termediate and the minor axes lie in PA~100°, and the 
line-tilt observed at PA=90° and 110° indicates that we 
are misaligned with both these axes. Thus, r^vpc—Ccspi 
along the apparent minor axis of the nebula, that is be¬ 
tween PA=90° and 110°. 

The cspl-zvpc connection for the inner and the outer 
shells is illustrated in Fig. 3, showing both r^vpc (at 
PA=90° and 110°) and Pexp vs IP. For reasons of ho¬ 
mogeneity each zvpc peak separation is normalized to the 
corresponding [O III] value. 

NGC 6818 is characterized by a large stratification of 
the radiation and of the kinematics, in agreement with the 
Wilson’s (1950) law. However, Pexp and r^vpc of [Ar IV] 
and [Ne III] considerably deviate from the sequence de¬ 
fined by the other ionic species. This behavior, also noticed 
in NGC 6565 (Paper IV), is typical of a PN powered by a 
high temperature central star (for details, see the above- 
mentioned reference and Sect. 9). Instead, the anomalous 
position of H I in Fig. 3 is the obvious consequence of the 
mono-electron atomic structure of hydrogen. 

The r^vpc vs IP relation being the same at all the ob¬ 
served PA (from Table 2), we can assess that in NGC 6818 
the expansion velocity is proportional to the distance from 
the central star through the relation: 

Vexp(fcms“^) = 3.5(±0.3) x r" (1) 

Weedman (1968) has derived Vexp{kms~^) =5.9 x 
(r" — 0.8), which is almost twice as steep as Eq. (1). The 
Weedman’s expansion law appears questionable, since it is 
based on the [O III] profile in a single, quite under-exposed 
Coude spectrum (R~60000) taken along the apparent ma¬ 


jor axis (PA~10°). Weedman noticed that A5007 A, al¬ 
though distorted by irregular motions, is un-tilted, and 
assumed “a priori” the following model: a prolate spheroid 
with a=17.3 arcsec and a/b=1.5, seen perpendicular to the 
major axis. 

Therefore in the following we will consider Eq. (1) as 
representative of the whole nebular kinematics. 

All this is synthesized in Fig. 4, showing the position- 
velocity (P-V) maps, i.e. the complete radial velocity field 
at the nine observed PA. They are relative to the systemic 
heliocentric velocity of the nebula, Vr0=-14.5(±l.O) km 
s“^, corresponding to Vlsr=-1-7(±1.0) km s“^, and are 
scaled according to Eq. (1). In other words, they repro¬ 
duce the tomographic maps in the nebular slices covered 
by the slit. We have selected He H, [O HI] and [N H] as 
representative of the high, medium and low ionization re¬ 
gions, respectively. 

These multicolor P~V maps highlight the large stratifi¬ 
cation of the radiation and the chaotic nebular structure. 
NGC 6818 is optically thin in most directions. The low 
excitation [N H] emission is associated to the presence 
of some inner and dense layer (see the moustaches, for 
example). The outermost, faint [N H] spikes detected in 
PA=10°, 90°, 150° and 170° belong to the external shell. 

The existence of a prominent hole to the North and of 
a smaller one to the South, as suggested by Rubin et al. 
(1998), is confirmed. The southern cavity is less evident 
in Fig. 4, being only grazed by our spectra at PA=170° 
and 10°. Further hollows are present (like at PA=70° and 
90°, in the approaching gas at West of the central star 
position). 

The radial ionization structure of the S-SE edge at 
PA=150°, corresponding to the region of the “cometary 
knot” visible in Fig. 1, is puzzling: [O HI] is lacking, 
whereas He H and [N H] are intense in the inner and the 
outer layers, respectively. 

In summary: although some spectral features of the 
northern and southern holes can be tentatively interpreted 
in terms of ionized gas moderately accelerated by some 
blowing-up agent (as suggested by Rubin et al. 1998), the 
linear expansion law here adopted represents a valid ap¬ 
proximation of the overall nebular kinematics. 

5. The physical conditions 

5.1. General considerations 

We can obtain the Te radial profile from the diagnostic 
line ratios of ions in p^ or p"' configurations (like [O HI] 
and [N H]), and the Ne radial distribution from both the 
diagnostics of ions in p^ configuration (like [S H]) and the 
absolute Ha flux. According to Paper IV, the zvpc, which 
is independent on the expansion velocity field, must be 
used. 

In the specific case of NGC 6818 the large Ha broad¬ 
ening (due to thermal motions, fine structure and expan¬ 
sion velocity gradient) prevents the accurate determina¬ 
tion of F(Ha)zvpc (and then of Ae(Ha)). Thus, in this 
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Table 2. Peak separation in the zvpc at the nine observed PA of NGC 6818 


Ion 2rzvpc (arcsec) 



PA=10° 

0 

O 

CO 

II 

0 

O 

II 

PA=70° 

0 

O 

II 

PA=110° 

PA=130° 

PA=150° 

PA=170 

[O I] 

23.5: 

23.0: 

22.0: 

—; 22.5: 

—; 22.5: 

—; 22.5: 

—; 22.0: 

—; 22.5: 

22.3: 

(S II] 

23.5 

22.4 

22.0 

15.6; 22.4 

—; 22.3: 

15.8; 22.5 

14.0; 22.0 

—; 22.4 

22.6 

(O II] 

23.1: 

22.8: 

21.7: 

— ; 22.2: 

—; 21.7: 

—; 22.5: 

—; 22.0: 

—; 22.0: 

21.8: 

HI 

19.2: 

16.7: 

16.4: 

14.0:; 21.0: 

12.8:; 19.5: 

13.0:; 21.0: 

13.1:; 20.5: 

14.7:; 21.0: 

20.2: 

[N II] 

23.5 

22.4 

22.0 

15.8; 22.4 

—; 22.2 

15.6; 22.4 

14.5; 22.4 

16.0; 22.4 

21.8 

(S III] 

21.3 

20.5 

18.4: 

15.3; — 

15.5:; — 

14.5:; — 

13.4; 21.3: 

15.0; 20.5: 

19.9: 

He I 

22.2 

20.5 

18.1 

15.7; 22.2 

14.5; 20.6 

14.5:; 21.0 

13.2:; 21.3 

—; 19.2: 

20.8 

(Ar HI] 

22.5 

19.8 

17.6 

15.4; 22.2 

14.6; 20.5 

14.7; 21.0 

13.7; 21.4 

15.4; 21.9 

21.1 

[0 HI] 

20.7 

20.0 

18.5 

15.3; 21.2 

15.0; 20.2 

15.0; 21.0 

13.8; 20.0 

15.9; 21.0 

20.4 

(Ar IV] 

17.9: 

15.3 

13.1: 

12.3; — 

11.0; — 

11.1; - 

12.0; — 

12.7; — 

16.4 

(Ne HI] 

21.7: 

20.9: 

20.0: 

15.3:; — 

15.2:; 21.0: 

14.5; 21.0 

14.3; 20.4 

—; 21.3: 

21.6 

He H 

19.0 

15.2 

12.3 

11.8; — 

11.0; — 

11.0; — 

12.1; - 

12.6; — 

16.6 

(Ar V] 

17.7: 

13.3: 

10.6: 

9.8; — 

8.0:; — 

8.5; — 

9.9; — 

10.9; — 

13.7 



' ' ' 1 ' ' ' 1 ' ' ' 1 ' ' ' 1 ' ' ' 

' ' ' ' 1 ' ' ' ' 

• PA=90 in - 


- 

- 

OPA=90 out- 

60 

X - 

m 

■ PA= 110 in— 


X 

• 

□ PA= 110 out 

eV] 

O 

- X X [ArlV] - 

[ArlV] • ^ [Nelll] - 

'—' 

X 

- 

• 

cu 

_ 

. 

_ 

1 _; 

X 

O 




20 

- 

- 

- 


^ X HI 

HI • ( 



X 

- 

■D 

0 

- X 

_1_1_1_1_1_1_1_1_1_1_1_1_1_1_1_1_1_1_l_ 

_1_1_1_1_1_1_1_1_1_ 

□0 

1 1 1 1 1 1 1 


80 60 40 20 0 0.5 1 1.5 

2Vexp [km/s] r(ion)/r([0III]) 


60 


40 


20 


0 


Fig. 3. Expansion velocity in the cspl (left) and peak separation in the zvpc for the inner and the outer shells of NGC 
6818 at PA=90° and 110° (right) vs IP. Each zvpc peak separation is normalized to the corresponding [O III] value. 
The ions in peculiar positions are marked. 


Section re[0 III], Te[N II] and Ne[S II] are derived from 
the corresponding line intensity ratios. Later on (Sect. 
6.3) we will illustrate the adopted escamotage providing 
F(Ha)zvpc and Ne(Jla) from the observed radial ioniza¬ 
tion structure and the assumption 0/H= constant across 
the nebula. 


5.2. Interstellar absorption 


First of all the observed line intensities must be corrected 
for interstellar absorption according to: 


log 


I(A)c 


1(A) 


obs 


= f\ c(H/3) 


( 2 ) 


where Ia is the interstellar extinction coefficient given by 
Seaton (1979). The logarithmic extinction at H/3, c(H/3), 





S. Benetti et al.: The Planetary Nebula NGC 6818 


7 


This figure is available in JPG format 


Fig. 4. Position-velocity maps at the nine observed PA for the high (He II, blue), medium ([O III], green) and low 
([N II], red) ionization regions of NGC 6818, scaled according to the relation Pexp(km s“^)= 3.5xr”. The orientation 
of these tomographic maps is the same of Fig. 2. 


is normally obtained by comparing the observed Balmer 
decrement (in particular Ha/H/3) to the intrinsic value 
given by Brocklehurst (1971) and Hummer & Storey 
(1987). The estimates of c(H/3) (from the Balmer ratio) 
reported in the literature for NGC 6818 span the range 
0.25 (Aller & Czyzak 1983, Liu & Danziger 1993) to 0.41 
(Collins et al. 1961). Moreover, Tylenda et al. (1992) and 
Condon et al. (1999) obtained c(H/3)=0.33 and 0.40, re¬ 
spectively, from the radio to H/3 fluxes. 

Thanks to the excellent spatial and spectral accura¬ 
cies achieved by the superposition technique used (±0.15 
arcsec and ±1.0 km s“^, respectively) we can extend the 
Ha/H/3 analysis to the whole spectral image, as recently 
introduced (Paper IV) in the study of NGC 6565, a com¬ 
pact, dust embedded PN exhibiting a complex c(H/J) pro- 
hle. The results are less dramatic for NGC 6818, since the 
blurred appearance of both Ho and H/3 limits the resolu¬ 
tion: besides some indications of a soft decline in the inner¬ 
most regions (likely caused by a local increase of Te), the 
spectral maps appear quite homogeneous at Ha/H/3=3.73 
(±0.06), corresponding to c(H/3)=0.37 (±0.03) (for the 
case B of Baker & Menzel 1938, re=12000 K and log Ne= 
3.00; Brocklehurst 1971, Aller 1984, Hummer & Storey 
1987). 

5.3. Te[0 III], Te[N II] and Ne]S II] 

First we obtain re[0 HI] from A5007 A/A4363 A, the line 
ratio being almost independent on Ne for Ve<10^ cm“^. 
Ne[S H] is then derived from A6717 A/A6731 A (using 
re[0 HI] to take into account the weak dependence of the 
ratio on Te). Last, Te[N H] comes from A6584 A/A5755 A 
(adopting Ne[S H] for its weak dependence on the electron 
density). 

Although the resulting Te[0 HI], Te[N H] and Ne[S H] 
profiles rapidly change with PA, as expected of the chaotic 


structure of NGC 6818, nevertheless there are some com¬ 
mon features. In order to highlight both the differences 
and the analogies we have selected two PA close to the 
apparent major axis (PA=10° and PA=30°) and two PA 
close to the apparent minor axis (PA=90° and PA=110°) 
as representative of the whole nebular phenomenology. 

The results are shown in Fig. 5; their main limitation 
is evident: due to the weakness of the [N H] auroral line 
and of the (S H] doublet, Te[N H] e Ne[S H] can be ob¬ 
tained only at the corresponding intensity peaks (in the 
best case). 

re[0 HI] presents a well-defined radial profile common 
to all the PA: it is >15000 K in the tenuous, innermost 
regions, it gradually decreases outward down to 12000- 
12500 K in the densest layers, and later it remains more 
or less constant (the last statement is weakened by the 
auroral line faintness). Such a radial trend is in quanti¬ 
tative agreement with the results by Rubin et al. (1998), 
based on A4363 A and A5007 A HST/WFPC2 imagery. 
The presence of a Te gradient across the nebula is also 
suggested by Hyung et al. (1999). Previous ground-based 
Te[0 HI] determinations are mean values for the bright¬ 
est (i.e. densest) regions and span the range 11400 K (de 
Freitas Pacheco et al. 1991) to 12770 K (Mathis et al. 
1998). 

In Fig. 5 Te[N H] refers to the intensity peaks of the 
low ionization regions and is systematically below Te[0 
HI], in agreement with both the previous reports for NGC 
6818 (9500 K, Hyung et al. 1999, to 11280 K, McKenna et 
al. 1996) and the general results for PNe (see Aller 1990, 
Gruenwald & Viegas 1995, Mathis et al. 1998). More Te 
values for NGC 6818 are: 11500 K (C HI], Kaler 1986), 
12700 K and 13800 K ([Ne V] and [O IV], respectively, 
Rowlands et al. 1989), 14700 K (Balmer discontinuity, Liu 
& Danziger 1993), 12500 K (C HI], Mathis et al. 1998), 
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Fig. 5. The diagnostics radial profile at PA=10° and 30° 
(close to the major axis of NGC 6818), and at PA=90° 
and 110° (close to the minor axis). Left ordinate scale: 
Te[0 III] (continuous line) and Te[N II] (squares). Right 
ordinate scale: 7Ve[S II] (circles). Te[N II] is lacking at 
PA=90° and 110°, because of the auroral line weakness. 


13000 K ([Cl IV], Hyung et al. 1999) and 19000 K ([O II], 
Keenan et al. 1999). 

In summary, the mean kinetic energy of the free elec¬ 
trons in NGC 6818 is quite large. This on the one hand 
explains both the Ha broadening and the strength of the 
forbidden lines (in particular, A5007 A of [O III]), on the 
other hand is indicative of a very hot central star (Th.> 
150000 K, also supported by the presence of high excita¬ 
tion emissions, up to [Ne V], IP=97.1 eV; see Table 1). 

Concerning the [S II] electron densities in the zvpc 
(Fig. 5), they are limited to the brightest parts of the 
external, low ionization layers and show peaks up to 
2000(±200) cm“^. Previously Hyung et al. (1999) ob¬ 
tained Ne[S H]~ 2000 cm“^ (they also report Ne[S 
H]~3000 cm“^ from the improved calculations for 3 equiv¬ 
alent p-electrons by Keenan et al. 1996). 

Besides the zvpc, we have extended the A6717 A/A6731 
A analysis to the prominent knots of the entire [S H] spec¬ 
tral images. The resulting iVe[S H] values span the range 
1500(±200) to 2800(±200) cm“^, the equatorial mous¬ 
taches being the densest regions of NGC 6818. 


In the next Section we will derive the local filling fac¬ 
tor, ei, by combining Ne[S H] and fVe(Ha)zvpc, since Ne[S 
Hjx e°-®~Ne(Ha)zvpc (Aller 1984, Osterbrock 1989). 


6. The radial ionization structure 

6.1. General considerations 

All the main ionic species are present in the echellograms, 
with the usual, dreadful handicap for hydrogen: due to the 
large Ha broadening, the deconvolution for instrumental 
resolution plus thermal motions plus fine structure ap¬ 
pears inadequate, and the detailed F(Ha)zvpc profile can¬ 
not be determined. Moreover: 

- the principal 0+ emissions (A3726 A and A3729 A) 
fall outside our spectral range. We have considered the 
much weaker A7319.87 A, the strongest line of the red 0+ 
quartet; 

- A3967 A of [Ne HI] being at the extreme blue edge of 
the frame, the quantitative Ne^"*" analysis appears quite 
uncertain. 

In absence of Ha we are forced to adopt A5007 A as ref¬ 
erence line, thus obtaining the radial ionization structure 
relative to 0++ according to: 


0 ++ 


f(A(A+^))zvpc 

F(A5007A)zvpc 


fiTe,Ne) 


(3) 


where a=b for the forbidden lines, and a=b-|-l for the re¬ 
combination ones. Note that the solution of the equations 
of statistical equilibrium requires the detailed knowledge 
of the plasma diagnostics, whereas Ne[S H] has been de¬ 
rived only at the intensity peaks of the outermost, low 
ionization layers. Luckily, in the Ne range here consid¬ 
ered f{Te,Ne) reduces to f(Te) for most ionic species. 
Moreover all the emissions, with the exception of A5876 A 
of He I and A4686 A of He H, are forbidden lines, whose 
emissivity is a direct function of Te. Thus, in first approx¬ 
imation we can put iVe=1.5xlO^ cm“^, constant across 
the nebula; this introduces uncertainties up to ±10% (for 
0+/0++ and Ar+^/0++), which do not modify the re¬ 
sults here obtained. 

The radial profiles of at the four selected PA of 
NGC 6818, shown in Fig. 6, contain a number of interest¬ 
ing features: 


- the weakness of the low excitation emissions in PA=10° 

northern sector, PA=90° both sectors and PA=110° 
both sectors indicates that the nebula is optically thin 
in these directions. It is almost thick in PA=10° south¬ 
ern sector and PA=30° both sectors; 

- as expected, He++/0++ decreases outward and He+/ 

0++ increases. They cross for He“''+/0'''+=He+/0+“'' 
~120, implying that Hetot/0++ ~240. Since in these 
internal regions Otot —1.30x0++ (see Eqs. (10) and 
(11) in Sect. 6.3) we derive He/0~185; 

- from similar considerations we infer N/Ocs 0.23, iVe/Oc:i 

0.15, S/O- 0.012 and Ar/O- 0.008. 
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Fig. 6. The radial ionization structure (relative to at PA=10° and 30° (close to the major axis of NGC 6818), 

and at PA=90° and 110° (close to the minor axis). The orientation is as in Fig. 5. 


These chemical abundances (relative to oxygen) must 
be compared with the corresponding values (relative to 
hydrogen) obtained from the conventional method, as il¬ 
lustrated in the next Section. 


6.2. Total chemical abundances 

According to the critical analysis by Alexander & Balick 
(1997) we consider the total line fluxes (i.e. integrated over 
the whole spatial profile and the expansion velocity held). 
The resulting ionic abundances must be multiplied for the 
corresponding ICFs, the correcting factors for the unob¬ 
served ionic stages. These were obtained both empirically 
(Barker 1983, 1986) and from interpolation of theoretical 


nebular models (Shields et al. 1981, Aller & Czyzak 1983, 
Aller 1984, Osterbrock 1989). 

The final mean chemical abundances of NGC 6818, 
presented in Table 3 (last column), are in reasonable agree¬ 
ment with the previous estimates reported in the literature 
(also listed in the Table), and in excellent agreement with 
the indications of Sect. 6.1. 


6.3. F(Ha)zvpc and Ne(Ha) 

As emphasized in Paper IV, the Ha flux distribution in 
the zvpc, F(Ha)zvpc 5 is the fundamental parameter linking 
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Table 3. Total chemical abundances (relative to hydrogen) 


Element 

Aller & Czyzak 
(1983) 

de Freitas Pacheco et al. 
(1991) 

Liu & Danziger 
(1993) 

Hyung et al. 
(1999) 

This paper 

He 

0.107 

0.126 

0.114 

0.105 

0.106(±0.003) 

C 

4.47xl0“'‘ 

- 

- 

8.0x10“"' 

- 

N 

1.41x10“^ 

1.1x10“"' 

- 

4.0x10“"' 

1.4(±0.2)xl0“"‘ 

0 

5.50x10“^ 

5.25x10“"' 

6.37x10“"' 

7.0x10“"' 

5.5(±0.5)xl0“"‘ 

Ne 

1.23xl0"'‘ 

- 

- 

1.0x10“"' 

9.0(±2.0)xl0“® 

Na 

3.09x10“® 

- 

- 

3.0x10“® 

- 

Mg 

- 

- 

- 

3.0x10“® 

- 

Si 

- 

- 

- 

9.0x10“® 

- 

S 

8.9x10“® 

1.2x10"^ 

- 

7.0x10“® 

6.2(±1.4)xl0“® 

Cl 

1.9x10“'^ 

- 

- 

3.0x10“'^ 

- 

Ar 

3.8x10“® 

- 

- 

4.0x10“® 

4.0(zb0.8)xl0"® 

K 

1.02x10“’^ 

- 

- 

2.0x10“'^ 

- 

Ca 

1.15x10“'^ 

- 

- 

1.5x10“'^ 

- 


Ne with both the spatial and the kinematical properties 
of the expanding plasma through the relation: 

re-O -47 ^eixrespixZ?^ 

where: 


- D is the nebular distance; 

- rcspi is the angular radius of the cspl (i.e. the nebular 

size in the radial direction); 

- e\ is the “local filling factor”, representing the fraction of 

the local volume actually filled by matter with density 
Ne. 


In order to recover F(Ha)zvpc we start assuming 0/H= 
5.5x10“'^, constant across the nebula. At each radial po¬ 
sition: 


H H° + H+ 


(5) 


It can be written in the form: 

O 0++ . 

- = -^x icf{0++) 


That is: 


O 

H 


F(A5007l)zvpc 
7^(77o;)zvpc 


X f{Te,Ne) x icf{0^'^) 


( 6 ) 

(7) 


Thus obtaining: 

7^(77o:)zvpc ~ 

= 1.9 X 103t(A5007A)zvpc x /(Te, Ne) x icf{0++) (8) 


Also in this case f{Te,Ne) essentially reduces to 
f{Te), with opposite trends for the emissivity of the for¬ 
bidden and the recombination line. 

Eq. (8) provides the Ha flux distribution in the 
zvpc (and fVe(Ha) through Eq. (4)) once the ioniza¬ 
tion correcting factor icf{0^'^) is known. The complex 


structure of NGC 6818 implies that ic/(0++) strongly 
changes across the nebula. In the innermost regions we 
have: 0°/0+'''<<l, 0+/0++<<l and H°/H+<<1. Thus 
ic/(0++) becomes: 


tc/(0++)in„er 


1 + 


0 ++ 


(9) 


Q+i 

According to Seaton (1968), ^q ++ — can be derived 
from the ionization structure of helium thanks to the 
closeness of the and He+ ionization potentials (54.9 
and 54.4 eV, respectively). To this end we have per¬ 
formed a number of photo-ionization simulations (see also 
Alexander and Balick 1997) obtaining the fairly good re¬ 
lation: 


0 ++ 


~ 0.30 X 


ffe++ 

He+ 


( 10 ) 


Therefore we will adopt: 


tc/(0++)in„er = 1 + 0.30 X —— (11) 

He~'' 

Concerning the outermost nebula, the contribution 
of J 2^=3 77'’"* can be neglected. Moreover in these re¬ 
gions we have 0°/0+<l (see Fig. 6), implying a low ef¬ 
ficiency of the charge-exchange reaction 0+ -I- 
-I- H+ (Williams 1973; Aller 1984; Osterbrock 1989), i.e. 
H°<<H+. Thus: 


*c/(0^~'’)outer — 1 + Q++ + 0++ 

All this is summarized in Fig. 7, showing the Ne ra¬ 
dial profile (at the four selected PA of NGC 6818) ob¬ 
tained from F(Ha)zvpc for some representative values of 
eix rcspi xD (in arcsec kpc), superimposed to A^e[S II] 
(taken from Fig. 5). 

Fig. 7 evidences the basic advantage of using 
F(Ha)zvpc in the determination of the electron density 
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Fig. 7. The iVe(HQ;)zvpc radial profile at the four selected PA of NGC 6818 for some representative values of ei x rcspi x D 
(dotted line=5 arcsec kpc, short-dashed line=10 arcsec kpc, long-dashed line=15 arcsec kpc), superimposed to A^e[S 
II] (thick line; same as Fig. 5). The orientation is as in Figs. 5 and 6. 


radial distribution: iVe(HQ;) extends all over the nebular 
image, whereas iVe[S II] is limited to the peaks of the low 
excitation regions. 

In detail: 

- at PA=10° (along the apparent major axis), A^e(Ha) 

presents a single, broad and asymmetric (i.e. steeper 
outwards) bell-shaped profile; 

- the double-peak structure is very subtle at PA=30° 

(close to the apparent major axis), the inner peaks 
being predominant; 

- at both PA=90° and 110° (close to the apparent mi¬ 

nor axis) the two-shell distribution clearly appears; 
the Ne(HQ!) top corresponds to the inner peaks at 
PA=110°, whereas at PA=90° the peaks are almost 
equivalent; 


- an outward, low density tail is present at all directions, 
extending up to about 15 arcsec from the star. 

Note in Fig. 7 the match between iVe(HQ;) and A^e[S II] 
for eixrcspixD~ 9.5(±1) arcsec kpc, implying that eixDci 
0.88 (±0.10) kpc (see Table 2 and Sect. 4). Although the 
assumption ei=l provides a lower limit of 0.9 kpc for the 
nebular distance, we can no longer delay a better quan¬ 
tification of this fundamental parameter. 

7. The nebular distance, mass and age 

The only individual distance reported in the literature 
(based on rough parameters for the star and the nebula) 
goes back to Gurzadyan (1970), who derived an indicative 
value of 2.2 kpc from the observed to expected size of the 
He^^ zone. 
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Table 4. Statistical distances of NGC 6818 


Author 

Distance (kpc) 

Method 

O’Dell (1962) 

1.7 

Shklovsky (1956), i.e. ionized mass=constant 

Cahn & Kaler (1971) 

1.82-2.38 

Shklovsky 

Cudworth (1974) 

1.95 

proper motions 

Milne & Aller (1975) 

2.07 

Shklovsky (radio) 

Cahn (1976) 

2.2 

Shklovsky 

Acker (1978) 

1.6 

published statistical distance scales re-calibrated with individual distances 

Maciel & Pottasch (1980) 

1.46 

ionized mass-radius relation 

Daub (1982) 

1.29 

ionized mass-radio surface brightness relation 

Phillips & Pottasch (1984) 

4.45 

observed vs predicted radio fluxes 

Maciel (1984) 

1.5 

ionized mass-radius relation 

Amnuel et al. (1984) 

0.95 

radio surface brightness-radius relation 

Kingsburgh et al. (1992) 

2.46 

Shklovsky 

Cahn et al. (1992) 

1.87 

ionized mass-surface brightness relation 

van de Steene & Zijlstra (1995) 

1.58 

radio continuum surface brightness temperature-radius relation 

Zhang (1995) 

1.68 

ionized mass-radius and radio cont. surface bright, temp.-radius relations 

Gorny et al. (1997) 

2.1 

theoretical evolutionary tracks vs observed nebular and stellar parameters 

Mal’kov (1997) 

2.5 

theoretical evolutionary age vs observed dynamical age 

Cazetta & Maciel (2000) 

>2.9 

Peimbert & Torres-Peimbert (1983) classicstellar massii^stellar gravity and 
temperatureh^stellar luminosity 

Bensby & Lundstrom (2001) 

2.23 

ionized mass-radius relation 

Phillips (2002) 

0.63 

radio surface brightness-radius relation 


As compensation the nebula is present in a good 20 
catalogues of statistical distances (listed in Table 4). They 
provide the following mean values: 

<D>(Shklovsky)~2.0(±0.4) kpc 
<D>(ionized mass-radius relation)~1.7(±0.4) kpc 
<D>(surface brightness-radius relation)~1.3(±0.4) kpc 
<D>(other methods)~2.8(±1.0) kpc 

To be noticed that the overall properties inferred in 
the previous Sections (in synthesis: NGC 6818 is an opti¬ 
cally thin, almost thin in some directions, quite “normal” 
PN) suggest that <D>(Shklovsky), <D>(ionized mass- 
radius relation) and <D>(surface brightness-radius rela¬ 
tion) are, at least grossly, reliable. 

In order to derive the dynamical parallax we have anal¬ 
ysed the first and second epoch (1998.30 and 2000.45, re¬ 
spectively) HST images of NGC 6818, searching for the 
angular expansion of the ionized gas; when combined with 
Eq. (1), it provides a reliable nebular distance (Reed et al. 
1999, Palen et al. 2002). Since the target is at the centre 
of the planetary camera (PC) chip in the 2000.45 images, 
whereas it is quite off-axis in the 1998.30 ones, the cor¬ 
rection for optical camera distortions was performed with 
the IRAF/STSDAS task “drizzle” (see Fruchter & Hook 
2002) using the Trauger coefficients. 

No apparent shift is obtained from the couples of [O 
III] and [N II] frames. 

According to Reed et al. (1999), we infer that the an¬ 
gular expansion of the gas is <1.2x10“^ arcsec, that is 
^ <5.5x10“^ arcsec yr“^. Since D(pc)=0.211[Fexp(km 
s-i)/^ (arcsec yr-i)], we derive D(NGC 6818)> 1300 pc. 


We have attempted a better quantification of D(NGC 
6818) through the interstellar absorption-distance rela¬ 
tion (Lutz 1973, Gathier et al. 1986, Saurer 1995) given by 
the field stars with accurate my, me, spectral type and lu¬ 
minosity class. They were selected using the SIMBAD fa¬ 
cilities of the CDS, Strasbourg Astronomical Observatory. 

The resulting Ay-distance law is presented in Fig. 8, 
where the small cluster of data at the right edge refers to 
the super-giant stars of NGC 6822 (Barnard’s galaxy), a 
nearby Ir galaxy projected at about 40 arcmin S-SE of 
our nebula. 

Fig. 8 indicates that: 

- the galactic absorption rapidly changes in direction 

(note the spread of the low-weight points). This is con¬ 
firmed by the inspection of both the ESO/SERC and 
the Palomar Schmidt surveys, showing a variable back¬ 
ground crossed by a series of extended, broad and faint 
emitting filaments; 

- close to the nebula direction (high-weight data). Ay 

quickly increases up to D~ 1.0-1.5 kpc, and later 
is constant (as expected of the large galactic lat¬ 
itude of the held, Jb Ilj ~18°). Note that for 
the super-giant stars of the nearby galaxy NGC 
6822 we have Ay(tot)=Ay(galactic)-l-Ay(NGC 6822), 
where Ay (NGC 6822) changes from star to star, and 
Ay(galactic)=constant. For Ay(NGC 6822)=0 we in¬ 
fer Ay(galactic)~ 0.80(±0.1), in excellent agreement 
with the literature reports (see Gallart et al. 1996, and 
Bianchi et al. 2001); 

- from Ay=2.18xc(H/3) (Acker 1978) we have Ay(NGC 

6818)=0.81(±0.06) and D(NGC 6818) >1.5 kpc. 
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Fig. 8. The interstellar absorption-distance relation for 
the field stars of NGC 6818. Three symbols are used, re¬ 
ferring to stars at different angular separation from the 
nebula. 

Although in Fig. 8 all the solutions for D(NGC 
6818)>1.5 kpc appear legitimate, both the large galac¬ 
tic latitude of the nebula (|b II| ~18°) and the low radial 
velocity relative to the Local Standard of Rest, Vlsr=- 
1.7 km s“^ (Sect. 4), decidedly favour the lowest values, 
i.e. 1.5 kpc< D(NGC 6818)<2.0 kpc (in agreement with 
the information previously obtained from the statistical 
distance scales). 

In the following we will adopt D(NGC 6818)= 1.7(1 q 2) 
kpc, corresponding to a distance from the galactic plane 
\z\ ~0.5 kpc. 

The local filling factor in the nebula is ei ~0.5. The 
ionized mass (obtained in different ways: from the H/3 
flux, the radio flux and the observed Ne distribution; Aller 
1984, Pottasch 1984, Osterbrock 1989 and Paper IV) re¬ 
sults to be Mion—0.13(±0.03) M©, and the kinematical 
age tkin oc r/Vexp ~ 2300(±300) yr. 

Mion is close to the total nebular mass, NGC 6818 
being an optically thin (almost thin in some direc¬ 
tions) PN, whereas tpin represents a lower limit to the 
actual age, tNGC esisi since the dynamical history of 
the gas is unknown. We can obtain a reliable esti¬ 
mate of tNGC 6818 by assuming a nebular ejection at 
Vexpsuperwind< RexpNGG 6818 , followed by a constant ac¬ 
celeration up to VexpNGC 6818 - Rexpsuperwind comes from 
the OH/IR sources, commonly regarded as the PNe pre¬ 
cursors (Habing 1996). For Vexpsuperwind =15(±5) km s“^ 
(Chengalur et al. 1993, David et al. 1993, and Sjouwerman 
et al. 1998) we derive Ingg 6818 = 35 00(±400) yr, i.e. our 
nebula is rather young. 

8. The central star 

As mentioned in Sect. 2, the exciting star of NGC 6818 
is the northern component of a visual binary system. 


The WFPC2 frames taken through the broad-band filter 
F555W (Ac= 5407 A, bandwidth= 1236 A) give my = 
17.06(±0.05) for the central star and my = 17.73(±0.05) 
for the southern red companion, where the unknown star 
colors are the main source of inaccuracy. Please note that 
for D=1.7 kpc we have (My)o(red companion) ~ -1-5.75, 
as expected of a late spectral type (G 8 to KO) Main 
Sequence star: a further, although weak, sign in favour 
of the adopted distance. 

According to Tylenda et al. (1993) and Feibelman 
(1994), the exciting star of NGC 6818 presents a weak 
emission line spectrum in both the UV and the optical re¬ 
gions. In order to derive the H I and He II Zanstra temper¬ 
atures, TzH I and TzHe II respectively, we have obtained 
the total H/3 and A4686 A nebular fluxes from the overall 
line profile at each PA, assuming a circular symmetry of 
the spectral image. We find log F(H/3)obs=-10.49(±0.03) 
mWxm“^ and F(A4686 A)/F(H/3)=0.65(±0.04), in excel¬ 
lent agreement with the values reported by Aller & Czyzak 
(1979), Kohoutek & Martin (1981), Webster (1983), Acker 
et al. (1991), and Hyung et al. (1999). 

The resulting Zanstra temperatures are log(TzH I)= 
5.20(±0.04) and log(TzHe H)=5.24(±0.04). We note in 
passing that TzH I~TzHe H; this is quite surprising for 
a high excitation, optically thin (almost thin in some di¬ 
rections) PN. Previous determinations (Harman & Seaton 
1966, Martin 1981, Gathier & Pottasch 1988 and Mal’kov 
1997) span the range 4.82 to 5.20 for log(TzH I), and 
5.03 to 5.30 for log(TzHe H). Moreover Pottasch & Preite- 
Martinez (1983) and Preite-Martinez et al. (1991) find log 
T*=4.93 and 5.02-5.05, respectively, from the energy bal¬ 
ance of the nebula. 

The stellar luminosities are log L*/L 0 (TzH I)=3.0 
(±0.1) and log L*/LQ(TzHe H)= 3.1(±0.1) (using the 
bolometric corrections by Schonberner 1981, and adopt¬ 
ing D=1.7 kpc). 

NGC 6818 results to be a rather young PN surrounding 
a very hot star of relatively low luminosity. This suggests 
that: 

- the central star is a hydrogen-burning post-AGB star 

(the helium-burning ones evolve too slowly), 

- the stellar mass, M*, is larger than the average value 

(~0.60 Mq) of the PNe nuclei (Blocker 1995 and ref¬ 
erences therein). 

The detailed comparison with the theoretical evolution¬ 
ary tracks by Schonberner (1981, 1983), Iben (1984), 
Wood & Faulkner (1986), Blocker & Schonberner (1990), 
Vassiliadis & Wood (1994) and Blocker (1995) confirms 
that no solutions are possible for the He-burning post- 
AGB stars (they are too slow), whereas the H-burning 
nuclei give M*=0.625-0.65 Mq. 

Moreover, the position of the central star of NGC 6818 
in the log L-log T diagram coincides with the beginning 
of the luminosity decline at the end of the shell nuclear 
burning. In this evolutionary phase all models present a 
large luminosity gradient with a deep minimum at log 
L*/Lq cs3.0, i.e. the value obtained for the central star of 
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NGC 6818. For example the 0.625 Mq H-burning nucleus 
by Blocker (1995) (showing an astonishing resemblance 
with our star) has log Lh</Lq= 3.387 at t=3457 yr, and 
log Lh</Lq= 2.846 at t=3585 yr, i.e. it drops in luminosity 
by a factor of 3.5 in 128 yr. 

And indeed, the central star of NGC 6818 is clearly 
visible in the earliest photographic reproduction of the 
nebula (Pease 1917; plate taken on July 12, 1912 with the 
Mount Wilson 60-inch reflector); the rough photometric 
sequence defined by the field stars gives mB = 15.2(±0.4). 
More stellar magnitude estimates reported in the litera¬ 
ture are: mpg=14 (Curtis 1918), mpg=15 (Berman 1937), 
mpv= 15.6 (Cudworth 1974), my— 14.9 (Martin 1981), 
mB> 15 (Shaw & Kaler 1985), mB=16.97 and mv=17.02 
(Gathier & Pottasch 1988), and mv=17.06 (this paper, 
referring to the epoch 1998-2000). 

In addition, the first spectroscopic observation of NGC 
6818 (star-l-nebula. Palmer 1903, referring to 1901) in¬ 
dicates “a pretty strong continuum spectrum extending 
out to A372, with the following bright lines: H/3, Hy, 
HJ...”. Palmer observed a large sample of nebulae with 
the Crossley reflector-|-plate spectrograph, reporting that 
“the faintest continuum spectrum photographed was that 
of a fifteenth magnitude star”. 

In spite of the photometric heterogeneity of these data, 
the suspicion of a “historical” decline for the nucleus of 
NGC 6818 appears legitimate (and robust). Also note that 
the magnitudes of the star just before the luminosity drop 
were: my ~mpv cild.h and mB — mpg ~14.0 (for D=1.7 
kpc and c(H/3)=0.37). 

In summary, all the evidence is that: 

- the quite massive central star is rapidly fading, 

- the thin-thick transition is hanging over the nebula. 

Within this scenario the following caveat is in order: 
because of the luminosity drop of the fast evolving star, 
the ionization and thermal structure of NGC 6818 are out 
of equilibrium. The recombination processes prevail, being 
faster in the high ionized species (Tylenda 1986, Marten 
& Szczerba 1997); this qualitatively explains the observed, 
unexpected closeness of the H 1 and He 11 Zanstra tem¬ 
peratures: the nebula is still optically thin (almost thin 
in some directions), but H/3 maintains a “longer memory” 
than A4686 A of the past stellar luminosity. Although the 
large Ne of the gas implies a short “time lag” of the main 
nebula (a few dozen years), nevertheless both TzHe 11 and, 
even more, TzH 1 are slightly over-valued. In other words: 
NGC 6818 is over-luminous with respect to the present 
UV stellar flux. 

The overall picture of our nebula shows a remarkable 
series of analogies with NGC 6565. According to Paper IV, 
NGC 6565 is a young (2300 yr), optically thick ellipsoid 
embedded in a large cocoon of neutral, dusty gas. It is in 
a deep recombination phase (started about 400 yr ago), 
caused by the luminosity drop of the massive powering 
star (M*~O.65M0), which has almost reached the white 
dwarf locus (log L*/L0C:i2.O, log T* ~5.08). 


All this suggests the evolutive contiguity of the two 
PNe, NGC 6818 being a bit older, but also a bit less “ma¬ 
ture” than NGC 6565, i.e. the central star of the latter 
is more massive. No doubt NGC 6818 will develop in the 
near future most of the observational characteristics al¬ 
ready present in NGC 6565 (and in the other recombining 
PNe contained in Tylenda 1986 and Corradi et al. 2000). 

We have also carried out a search for other PNe in the 
peculiar evolutionary phase of NGC 6818 (e.g. at the very 
beginning of the thin-thick transition) adopting the follow¬ 
ing criteria: high temperature and quite low luminosity of 
the star, different [O HI] and [N H] nebular morphology, 
high excitation class and relatively weak low ionization 
emissions. Among the dozen or so candidates we identi¬ 
fied, the most promising are NGC 6326, NGC 6884, NGC 
7354 and IC 4663. 

9. The photo-ionization model 

9.1. General approach 

In this Section the photo-ionization code CLOUDY 
(Ferland et al. 1998) is applied to a PN having the same 
distance, matter distribution and chemical composition of 
NGC 6818. 

We immediately point out that CLOUDY is a “steady- 
state” model (a constant UV flux hits the nebula), whereas 
we argued that the ionization and thermal structure of 
NGC 6818 are out of equilibrium, because of the fast lu¬ 
minosity drop of the star. A similar, but less dramatic 
situation occurs in the recombining PN NGC 6565, which 
is in quasi-equilibrium, since the stellar evolution slows 
down and the luminosity gradient rapidly decreases while 
approaching the white dwarf domain (Paper IV). 

In absence of an “evolving” photo-ionization code, 
which takes into account the gas reactions to a fast 
changing UV flux, at the moment we apply CLOUDY 
to the present nebula, (NGC 6818)(to), as well to (NGC 
6818)(to-100yr) and (NGC 6818)(to-l-100yr). Both the 
past and the future “snapshots” come from a homolo¬ 
gous expansion (i.e. roct and A^(H)(xt“^), being the dy¬ 
namical time to= 2300 yr and At<<to. Moreover, the 
evolutionary track of the 0.625 M© H-burning post-AGB 
star by Blocker (1995) was adopted for the stellar pa¬ 
rameters. Last, we have selected the quite thin East sec¬ 
tor at PA=90° as “mean” density radial profile of (NGC 
6818)(to). 

The input data are contained in Table 5. The results, 
shown in Fig. 9, can be synthesized as follows: 

- at (to-lOOyr) (upper row) the intense UV flux of the 

luminous and hot star completely ionizes the nebula, 
which appears as a very high excitation object at large 
Te, optically thin in all the directions; 

- at (to) (central row) the stellar luminosity is hardly suf¬ 

ficient for ionizing the gas, and the high excitation neb¬ 
ula is almost thick in some directions, where the low 
ionization emissions emerge; 
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Fig. 9. “Steady” photo-ionization model of NGC 6818, East sector in PA=90°, at three epochs: (to-lOOyr) (upper row), 
to (central row), and (to-l-100yr) (bottom row), to=2300 yr being the dynamical nebular age. Left panels: radial profiles 
of the absolute line fluxes. Right panels: the diagnostics; left ordinate scale: Te (dotted line); right ordinate scale: Ne 
(thick continuous line), iV(Htot) (dotted-dashed thin line), (long-dashed thin line) and A^(H°) (short-dashed 

thin line). 


- at (to-l-100yr) (bottom row) the UV flux is far inade¬ 
quate: the innermost plasma maintains a high degree 
of excitation, but the external layers are thick, with 
prominent [O I], [O II], [N II] and ]S II] lines. Note that 
the outermost parts of the “steady” photo-ionization 
model are neutral, whereas are recombining in a “true” 
evolving nebula (they create a faint halo embedding 
the main object). 

Although these results call for caution due to the 
photo-ionization model limitations, nevertheless they sup¬ 


port the evolutive scenario of the previous Sections: NGC 
6818 is a PN at the very beginning of the recombination 
phase. 

A direct confirmation is expected of the detailed analy¬ 
sis of the equatorial moustaches, representing the densest 
and brightest (i.e. thickest) regions of NGC 6818 (Sect. 
5.3). To this end we have first de-projected the spectral 
images through Eq. (1), and then repeated the complete 
procedure already applied to the zvpc. 
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Table 5. The CLOUDY photo-ionization code: input parameters of the model nebula 


Radial density profile 


right panels of Fig. 9 (cfr. Sects. 6.3 and 9) 


Chemical abundances: 

C, Na, Mg, Si, Cl, K, Ca Hyung et al. (1999) 
He, N, O, Ne, S, Ar this paper 

other elements 


Dust 


PN (CLOUDY default) 
PN (CLOUDY default) 


Filling factor 


Exciting star 


0.5 


at time (to-lOOyr): blackbody with T*=170000 K and log L*/L 0 = 3.45 (*) 
at time (to): blackbody with T* = 160000 K and log L*/L 0 = 3.10 (see Sect. 8 ) 
at time (to-|-100yr): blackbody with T* = 140000 K and log L*/L 0 = 2.60 (*) 


(*) according to Blocker (1995) 


9.2. Equatorial moustaches 

Let’s consider the southern, approaching moustache in 
PA=10°. The de-projected radial prohles of the main 
emissions, the diagnostics and the ionic abundances 
are shown in Fig. 10. 

We are particularly interested to A6300 A of [O I], 
marking the neutral (or almost neutral) nebular layers, 
where the efficiency of the resonant charge-exchange re¬ 
action O’*" -I- -I- H+ largely increases (Williams 

1973). Note that on the one hand the [O I] line is very sen¬ 
sitive to the physical conditions, on the other hand no pre¬ 
cise information is yet available for the electron tempera¬ 
ture of the external, neutral (or almost neutral) gas dom¬ 
inated by the recombination and cooling processes. We 
arbitrarily adopt re(almost neutral gas)=8000 K, which 
is lower than Te(ionized gas), but not enough low to com¬ 
promise the [O I] emissivity (a choice supported by the 
large kinetic energy of the free electrons in NGC 6818; see 
Sect. 5.3). 

Moreover, referring to the ionization correcting fac¬ 
tor fc/(0++)outer of Scct. 6.3, it includcs H°, whose con¬ 
tribution comes from the equilibrium condition 0“'"/0 ° 
~0.82x(H+/HO), valid for 0°>0+, i.e. in the nebular 
regions affected by the O’*" -I- -I- H+ reaction 

(Williams 1973, Stancil et al. 1999). 

The observational results of Fig. 10 must be compared 
with the “steady” photo-ionization model of the mous¬ 
tache presented in Fig. 11 (the stellar parameters obvi¬ 
ously refer to time to of Table 5). 

Besides the general features common to all the direc¬ 
tions (i.e. the decreasing re[0 III] radial profile, re[N 
II] <Te[0 III], He++/0++=He+/0++ -120, IVe(Ha)= 
Ne[S II] for (ei x rcspixD) — 9.5 arcsec kpc, and so on), the 
moustache of NGC 6818 (Fig. 10) exhibits a remarkable 
peculiarity: the external layers are partially neutral, and 
the outward ionization decline is very smooth, as expected 
of a recombining region. On the contrary, the “steady” 


photo-ionization model in Fig. 11 predicts an abrupt ion¬ 
ization fall (to be noticed: the photo-ionization model at 
(to-lOOyr), not shown here, indicates that a century ago 
the moustache was optically thin to the UV stellar flux). 
The Ne depletion rate for recombination is given by: 

dNe/dt = -OB xNex N{H+) (13) 

where ob is the recombination coefficient (Storey & 
Hummer 1995). Integrating Eq. (13) and assuming 
Ne=1.15xN(H+), we obtain: 


Ne{to -t) - Ne{to) 

Ob X Ne{to — t) x Ne{to) 


(14) 


which provides the time t elapsed from the beginning of 
the recombination once are known Ne(to), the present 
electron density, and Ne{to — t), the electron density at 
the start of the process. 

The application of Eq. (14) to the external parts of the 
“true” moustache (Fig. 10) furnishes t=30-60 yr (t=60- 
120 yr and 20-40 yr for Te(almost neutral gas)=5000 K 
and 12000 K, respectively). In spite of the heavy assump¬ 
tions, this agrees with all the previous evidences suggest¬ 
ing that the recombination phase has just begun in NGG 
6818. 


10. The 3-D morpho-kinematical structure 

10.1. General 

The reconstruction of the gas distribution in the nebular 
slices covered by the slit was introduced in Papers I and 

11 . In the case of NGG 6818 we have selected A4686 A of 
He H, A5007 A of [O HI] and A6584 A of [N H] as rep¬ 
resentative of the high, mean and low ionization regions, 
respectively. Note that Ha, the marker of the whole ion¬ 
ized gas distribution, cannot be utilized because of the 
blurred appearance. 
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Fig. 11. “Steady” photo-ionization model for the south¬ 
ern, approaching moustache of NGC 6818 in PA=10°. Top 
panel: absolute radial flux distribution (erg cm“^ s“^) in 
the different emissions; same symbols as Figs. 6, 9 and 
10. Bottom panel: diagnostics radial profile; left ordinate 
scale: dotted line= Te; right ordinate scale: thick contin¬ 
uous line= Ne, long-dashed line= short-dashed 

line= iV(H°), dotted-dashed line= IV(Htot)- 


Fig. 10. The observed radial properties of the southern, 
approaching moustache of NGC 6818 in PA=10°. Top 
panel: the absolute flux (in erg cm“^ s“^) of the main 
emissions; same symbols as Figs. 6 and 9. Middle panel: 
the diagnostics; left ordinate scale: dotted line= Te[0 
III], short-dashed line= re[N II]; right ordinate scale: 
thick continuous line= A^e[S II], long-dashed line= 7V(H+) 
for (eiXr(.spixD)= 9.5 arcsec kpc, dotted-dashed line= 
A^(IItot). Bottom panel: the ionic abundances (same 
symbols as the top panel). 


The spectral images of the forbidden lines have been 
de-convolved for seeing, spectral resolution and thermal 
motions, while also fine structure has been taken into ac¬ 
count for A4686 A of He II. They all were de-projected 
through Eq. (I), and assembled by means of the 3-D ren¬ 
dering procedure described in Paper III. 

In order to reproduce the spatial structure on the pa¬ 
per we adopt the usual method: a series of opaque re¬ 
constructions of the nebula seen from different directions, 
separated by 15°. Each couple forms a stereoscopic pair 
providing a 3-D view of NGC 6818. 

The novelties are represented by the multicolor projec¬ 
tion and the movies. 


10.2. Opaque reconstruction 

For reasons of space we only present the opaque recon¬ 
struction in He H, [O HI] and [N H] for a rotation of 180° 
around the East-West axis (close to the minor axis). This 
is shown in Figs. 12 to 16, where the upper-right panel 
corresponds to the nebula seen from the Earth (West is up 
and North to the left, to allow the reader the stereo-view). 

The high excitation layers of NGC 6818 (Fig. 12) form 
an inhomogeneous shell seen almost equatorial on, open- 
ended at North and South; it merges in a closed ellipsoid 
at lower A4686 A fluxes. 

The opaque reconstruction in [O HI] is given for 
two absolute flux cuts: log E(A5007 A)=-17.32 erg s“^ 
cm“^ (Fig- 13), and -17.75 erg s“^ cm“^ (Fig. 14). Since 
O/H=5.5 x 10-4~O++/H+ and Ne=1.15xN(il+), they 
correspond to A^e~1500 cm“^ (Fig. 13), and A^e~900 
cm-3 (Fig. 14) (for Te=12000 K and ei=l). The [O HI] 
high cut represents the densest regions of the inner shell, 
mainly constituted by the equatorial moustaches and by 
an extended cup in the southern part, whereas the exter¬ 
nal, inhomogeneous shell appears in the low cut frames, 
characterized by the large hole at North (along the major 
axis). 

The brightest and the faintest [N H] regions refer to 
A6584 A at the cuts log E(A6584 A)=-18.20 erg s“^ cm“^ 
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Fig. 12. The structure of NGC 6818 for a rotation around the East-West axis centered on the exciting star. Opaque 
reconstruction at A4686 A of He II, seen from 13 directions separated by 15°. The line of view is given by {Oj'ip), 
where 9 is the zenith angle and ip the azimuthal angle. Each horizontal couple represents a “direct” stereoscopic pair, 
and the whole figure provides 12 3-D views of the nebula in as many directions, covering a straight angle. Following 
Paper III: to obtain the three-dimensional vision, look at a distant object and slowly insert the figure in the field of 
view, always maintaining your eyes parallel. Alternatively, you can use the two small dots in the upper part of the 
figure as follows: approach the page till the two dots merge (they appear out of focus); then recede very slowly, always 
maintaining the two dots superimposed, till the image appears in focus. The upper-right image is the rebuilt-nebula 
seen from the Earth (West is up and North to the left, to allow the reader the stereo-view). 


(Fig. 15), and -18.80 erg s“^ cm“^ (Fig. 16). The [N II] 
(high cut) distribution mimics the corresponding [O III] 
one, whereas at lower fluxes the outermost, knotty struc¬ 
ture emerges. Note the general weakness of the low ion¬ 
ization layers for -135° < ip < —45°. 

These stereo-reconstructions confirm that NGC 6818 
consists of a double shell projected almost equatorial 
on; the external one is spherical (r~0.090 pc), faint and 
patchy. It circumscribes a dense and inhomogeneous tri- 
axial ellipsoid (a/2~0.077 pc, a/b~1.25, b/c~1.15) char¬ 
acterized by a large hole along the major axis and a couple 
of thick equatorial regions (the moustaches). 

10.3. Multicolor projection 

A representative sample of the different morphologies as¬ 
sumed by NGC 6818 when changing the line of view 
is given in Figs. 17 and 18. They show the assembled, 
multicolor projection at high (He H, blue), medium ([O 
HI], green) and low ([N H], red) ionization for a rota¬ 
tion around the N-S axis (close to the major axis; Fig. 
17), and around the E-W axis (close to the minor axis; 
Fig. 18). Also in this case the upper-right panel corre¬ 
sponds to the nebula seen from the Earth (North is up and 
East to the left), to be compared with Fig. 1, and with 
the HST multicolor image by A. Hajian fc Y. Terzian at 
http://ad.usno. navy, mil/pne/gallery.html. 

These “almost” true color reproductions highlight the 
variety of looks exhibited by the Little Gem: roundish to 
elliptical, to quasi-bipolar (according to the current mor¬ 


phological classifications; Greig 1972, Stanghellini et al. 
1993, Corradi & Schwarz 1995, Gorny et al. 1997). NGC 
6818 resembles NGC 6153, Hu 1-1 and NGC 4071 when 
seen from (0,60), K 3-57 and M 2-51 from (0,90), A 70 
from (90,0), IC 4663 from (0,30) and NGC 7354 from 
(150,0) (see the imagery catalogues of PNe by Acker et 
al. 1992, Schwarz et al. 1992, Manchado et al. 1996 and 
Gorny et al. 1999). 

10.4. Movies 

A number of limitations are implicit when showing the 3- 
D structure on the paper. They are both objective (space, 
choice of the rotation axes, the nebular parameters, the 
cuts etc.), and subjective (difficulty in the stereo view). 
In order to overcome all these handicaps, we have decided 
to introduce a series of movies as integral and functional 
part of each PN analysis. 

Such an “in fieri” film library, providing the mul¬ 
ticolor projection, the opaque recovery in different 
ions and at various cuts, some slices and radial pro- 
files etc, (suggestions are welcome) , can be found at 
http://yyeb.pd.astro.it/sabbadin] (where the acronym 
“sabbadin” stands for “stratigraphy and best boundary 
analytic determination in nebulae”, sic!). 

We are also exploring new graphical solutions, render¬ 
ing at best the (present) spatial reconstruction and the 
(forthcoming) spatio-temporal one. In particular, we are 
sounding the wide potential of the virtual reality. 
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Fig. 13. Same as Fig. 12, but for [O III] at the flux cut log E(A5007 A)=-17.32 erg s ^ cm ^ {Ne^ 1500 cm ^ for 
Te=12000 K and ei=l). 


This figure is available in JPG format 


Fig. 14. Same as Fig. 13, but for the [O III] flux cut log E(A5007 A)=-17.75 erg s ^ cm ^ (A^e~ 900 cm ^ for 
re=12000 K and ei=l). 


11. Discussion and concluding remarks 

In the previous Sections of this paper dedicated to the 
Little Gem we have (partially) extracted and interpreted 
the huge amount of physical information contained in the 
ESO NTT+EMMI echellograms by means of the 3-D pro¬ 
cedure developed in Papers I to IV for all types of ex¬ 
tended, expanding nebulae. 

NGC 6818 results to be a young (3500 yr), optically 
thin (quasi-thin in some directions) double-shell at a dis¬ 
tance of 1.7 kpc, projected almost equatorial-on: a ten¬ 
uous, patchy spheroid (r~0.090 pc) encircles an inner, 
dense, inhomogeneous ellipsoid (a/2~0.077 pc, a/b~1.25, 
b/c~1.15) empty along the major axis and optically thick 
in a pair of equatorial moustaches. 

We recall that outer attached shells have been found 
in both the ID and 2D hydrodynamical models (see e. g. 


Mellema 1995, Corradi et al. 2000 and Villaver et al. 2002), 
due to a D-type ionization front in the early evolution 
of the PN. However, all these models fail to reproduce 
the smooth, innermost radial density profile observed in 
NGC 6818 (Fig. 7); they predict an empty region (the hot 
bubble) and a quick density rise (the gas compression is 
provided by the thermal energy of the hot bubble formed 
by the adiabatic shock at the interaction region between 
the high velocity stellar wind and the material ejected 
during the superwind phase). 

The central star of NGC 6818 is a visual binary: a 
faint, red companion appears at 0.09 arcsec in PA=190°, 
corresponding to a separation >150 AU. For an orbit of 
low eccentricity the Kepler’s third law provides a period 
>1500 yr. Note that (by chance?) the two stellar compo¬ 
nents appear aligned with the major axis of the nebula. 
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Fig. 15. Same as Figs. 12 to 14, but for [N II] at the flux cut log E(A6584 X) 
ionization layers). 


18.20 erg s ^ cm ^ 


(bright low 


This figure is available in JPG format 


Fig. 16. Same as Fig. 15, but for the [N II] flux cut log E(A6584 A)=-18.80 erg s ^ cm ^ (faint low ionization layers). 


Despite some pioneering studies, the physical effects 
produced by a wide binary system on the PN ejection and 
shaping are still poorly known. Following Sober (1994), 
an orbital period comparable or longer than the mass-loss 
episode generating the nebula causes a density enhance¬ 
ment in the equatorial plane and/or spiral structures. 
Soker (2001) suggests that in wide binary systems (final 
orbital periods in the range 40 to 10'^ yr) an outer, spheri¬ 
cal structure is formed by the early AGB wind. Toward the 
end of the AGB phase, the increased mass-loss rate creates 
an accretion disk around the companion. If this blows jets 
or a collimated fast wind, two lobes appear in the inner 
nebula (a multi-lobed structure in the case of a processing 
accretion disk). Always according to Soker (2001), a frac¬ 
tion of 5% to 20% of all PNe originate in such wide binary 
systems. Curiously, the same author (Soker 1997) includes 
(with a high degree of confidence) NGC 6818 among the 
PNe resulting from the common envelope evolution with a 


sub-stellar companion (planet(s) and/or a brown dwarf). 
It is evident that the argument deserves further attention. 

The Little Gem is in a peculiar evolutionary phase, 
i.e. at the very beginning of the recombination process. 
This is caused by the luminosity decline of the 0.625 M© 
central star (log T* ~5.22 K; log L^/Lq ~3.1), which has 
recently exhausted the nuclear shell burning and is rapidly 
moving towards the white dwarf region. The stellar drop 
being fated to continue, NGC 6818 will become thicker 
and thicker, and the amount of neutral, dusty gas in the 
outermost layers will increase with time. The ionization 
front will re-grow only in some hundreds years, when the 
gas dilution due to the expansion will overcome the slower 
and slower luminosity decline. 

Concerning the observational analogies between NGC 
6818 (this paper) and NGC 6565 (Paper IV), in Sect. 8 
we have pointed out the probable evolutive contiguity of 
the two nebulae. We stress here a more facet of the af- 
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Fig. 17. Multicolor appearance of NGC 6818 (blue=He II, green=[0 III], red=[N II]) for a rotation around the N-S 
axis centered on the exciting star. The right panel, (0,0), corresponds to the re-built nebula seen from the Earth (North 
is up and East to the left). Same scale as Figs. 12 to 16. Recall that projection(0, ^)=projection(0 ± 180°, ^/>± 180°). 


This figure is available in JPG format 


Fig. 18. Same as Fig. 17, but for a rotation around the E-W axis. 


fair, e.g. the importance of the temporal factor: thanks 
to the excellent spatial and spectral resolutions achieved 
by the 3-D analysis applied to high quality spectra, we 
can no longer regard a PN as a static, uniform and un¬ 
changeable object; it is a dynamical, inhomogeneous and 
evolving plasma. 

Ironically, this quite reverses the gap between theory 
and practice outlined in the Introduction: now the existing 
“steady” photo-ionization models appear inadequate to 
interprete the observational data. A 3-D “evolving” code 
is highly desired, providing for the gas reactions to the 
changing UV flux of the central star. 

In summary: we have inferred a self-consistent picture 
of the Little Gem by means of ESO NTT-I-EMMI echel- 
lograms. Deeper observations at even higher spatial and 
spectral resolutions will disentangle the still unresolved 
problems, like the accurate Te[N II] and A^e[S II] radial 


profiles (Sect. 5), the intriguing ionization structure of 
the cometary knot in PA=I50° (Sect. 4) and the possi¬ 
ble blowing of the gas in the Northern and Southern holes 
(Sect. 4). Moreover, a gradual change of the [N II]/[O III] 
morphology is expected in the future HST imagery, due to 
the peculiar evolutionary phase of the nebula. Concerning 
the exciting star, a painstaking search in the world-wide 
archives (both spectroscopic and photometric), and new, 
deep, UV to IR spectra of the stellar system (hot central 
star -I- cold companion) are needed. 

At last, the “vexata quaestio”: which are the mecha¬ 
nisms and the physical processes ejecting and shaping a 
PN like NGC 6818? In our opinion the question appears 
premature, and the answer is beyond the aims of the pa¬ 
per, given the “forest” of proposed models (see Icke et al. 
1992, Mellema 1997, Dwarkadas & Balick 1998, Garcia- 
Segura et al. 1999, Frank 1999, Blackman et al. 2001, 
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Soker & Rappaport 2001, Balick & Frank 2002), and the 
“desert” of carefully studied true nebulae. We are confi¬ 
dent that new, reliable and deep insights on each object 
and on the whole class will come out of the comparative 
analysis of a representative sample of PNe and proto-PNe. 

This is the final goal of our survey carried out with 
ESO NTT-bEMMI and TNG-bSARG. Indeed, the superb 
quality of these echellograms constitutes a powerful tool 
for unveiling the evolutional secrets of the PNe (as well 
for masking the cultural gaps of the authors). 
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